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^_^ Abstract. We consider the Rcnyi entropies Sn in one-dimensional massive 

^ integrable models diagonalizable by means of corner transfer matrices (as 

Heisenberg and Ising spin chains). By means of explicit examples and using the 
relation of corner transfer matrix with the Virasoro algebra, we show that close 
to a conformal invariant critical point, when the correlation length ^ is finite but 
^H la^rge, the corrections to the scaling are of the unusual form ^—^/"^ with x the 

, 1^ dimension of a relevant operator in the conformal theory. This is reminiscent of 

^ the results for gapless chains and should be valid for any massive one-dimensional 

^— * model close to a conformal critical point. 
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1. Introduction 

In the last decade there has been an increasing interest in quantifying the amount 
of entanglement that is present in the ground state of an extended quantum system 
[T]. A measure of the bipartite entanglement is given by the so-called Renyi entropies, 
defined as follows. Let |^) be the ground state of an extended quantum mechanical 
system and p — |^)(^| its density matrix. One then divides the Hilbert space into a 
part A and its complement B and considers the reduced density matrix p^ = Tr^ p of 
subsystem A. Finally, the Renyi entropies are given by 

Sn^^lnTrpl. (1) 

The particular case n = 1 of (fTl) is the von Neumann entropy 5*1 and it is usually 
called simply entanglement entropy. However, the knowledge of Sn for different n 
characterizes the full spectrum of non-zero eigenvalues of p^ (see e.g. |2]) and provides 
significantly more information than Si . 

The case that has been most widely studied is a critical one-dimensional system 
whose continuum limit is described by a conformal field theory (CFT) of central charge 
c. When A is an interval of length £ embedded in an infinite system, the asymptotic 
large-^ behavior of the Renyi entropies is given by [2 HI [S] 

^nW^^(^l + ^)ln^ + c;, (2) 

where c'„ is a non-universal constant. This behavior has been verified analytically and 
numerically in numerous models. Only more recently it has been realized that such 
asymptotic behavior can be difficult to observe because of the presence of large and 
unusual corrections to the scaling behaving like [SI [71 15] 

Sn{e)c^^(l + ^yn£ + c; + 6„f-2./„ ^ (3) 

where x is the scaling dimension of a relevant operator (i.e. x < 2) and 6„ another non- 
universal constant. These power-law corrections decay slowly even for moderate values 
of n and completely obscure the asymptotic result for large n. In field theory language, 
this unusual scaling originates from the conical singularities (at the boundary between 
A and B) present in the n-sheeted Riemann surface describing Trp^ [7]. 

However, universal scaling is not a prerogative of the gapless models. A nearby 
critical point influences a part of the parameter space called 'critical region' in which 
the correlation length ^ (inverse gap) is large but finite. Simple scaling arguments 
suggest that when an infinite system is divided in two semi-infinite halves, the 
entanglement entropy should scale as 

^n^^(l + ^)lne + C;, (4) 

where C'^ is yet another non-universal constant. (One could also fix the normalization 
of the correlation length in such a way that C^ = c^/2 as done in some practical 
instances |9j.) This formula has indeed been corroborated by a general field-theory 
argument that parallels the c-theorem [4] and it has been verified for many integrable 
models. In this massive case, the spatial division in A and B is not as important 
as in the massless case. Indeed as long as the correlation length is smaller than all 
the separations (e.g. £ above), the Renyi entropies of many disjoint blocks are given 
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by Eq. (HI) multiplied by the number of boundary points between A and B. This 
is nothing but the one-dimensional area law. The question that naturally arises is 
whether unusual corrections to the scaling in ^ would be also present for the massive 
case. By simple scaling hypothesis, one could argue that corrections due to a finite 
correlation length in Eq. (Is]) should be functions of S,/i and then changing the limit 
from £ -^ £^ to i ^ £^, one argues 

^„=^^(l + ^)ln? + C; + S„r"/". (5) 

(Notice that the exponent of the corrections is half the one in (Is]) because of the 
presence of a single conical singularity, analogously to the case of critical systems with 
boundaries [7j). However, because of the very unusual form of these corrections, one 
could doubt whether such a simple scaling argument gives the correct answerjf] Here 
we provide the analytical evidence that in some integrable models in which Sn can be 
obtained through the Baxter corner transfer matrix [TU] this scaling is indeed correct 
and that further corrections are of the form ^-^^/^ and ^-^^ with k integer. 

2. Corrections to the scaling of Renyi entropies using corner transfer 
matrix 

When dealing with the geometric bipartition considered in this paper (i.e. two 
semi- infinite half lines) the corner transfer matrix (CTM)- that is a classical tool of 
statistical mechanics [10]- helps considerably in the derivation of the reduced density 
matrix p^ and hence for the characterization of the entanglement entropies [H [TT] . 

In two-dimensional statistical models the CTM A is the transfer matrix that 
connects an horizontal row (let say the line a; < 0) to a vertical one (let say y < Q). 
Rigorously speaking, each element of A is the partition function of the left-lower 
quadrant when the spins on the negative x and y axes are fixed to given values. If, 
following Baxter [TU], we choose the lattice in a clever way (i.e. rotated by 7r/4 with 
respect to the axes), the four corner transfer matrices are all equivalent (under some 
symmetry requirements on the model) and the full partition function is just TrA^. 
It is also evident that the reduced density matrix of the quantum problem whose 
time-discretized version is the classical model under consideration (see e.g [TU[T3]) is 

^-^ = ^' ^ ^ (6) 

where the denominator is fixed by the normalization Trp^ — 1. 

For several models, it is possible to obtain this reduced density matrix exactly in 
the infinite-lattice limit.. It always assumes the form [131 [TT] 

p — -^CTM 

where -ffcTM is an effective Hamiltonian, which, for many interesting integrable 
massive models, may be rewritten in terms of free-fermionic operators (see |13j and 
references therein) 



3=0 



HcTM = V CjUj . (8) 



I This scaling analysis is corroborated by recent results for gaploss systems in a trapping potential 
[9|. In this case, the actual system is gapless, but a length scale is generated by the trapping potential 
and the resulting scaling is exactly given by Eq. (|5|l where J is replaced by the trap-scale. 
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where Uj — clci are fermion occupation numbers and have eigenvalues and 1. These 
properties are very general features of a large class of massive integrable spin-chains. 
The specificity of the model comes from the functional form of the "single-particle 
levels" Ej. In this basis p^ is diagonal and so the Renyi entropies are easily written 
for any value of n, not only integer. 

One should mention that, for a given lattice Hamiltonian, there are several ways 
to construct a CTM. The one used above is particularly convenient because then the 
RDM is just its fourth power. In general, pA will be given by the product 'ABCD' 
of four different CTMs which also not need to be symmetric (although the product 
is). CTMs can be defined also for non- integrable models, but then one does not have 
closed and simple formulae for the eigenvalues and has to work numerically, see e.g. 

Ref. ng. 

2.1. Non-critical XXZ chain 

We start our analysis from the anisotropic Heisenberg model with Hamiltonian 

j 
where cr's are the Pauli matrices. This model is gapless for |A| < 1 and gapped 
for |A| > 1 with a conformal point at A — 1. We study the antiferromagnetic 
gapped regime with A > 1. The XXZ model is usually solved by Bethe ansatz 
that is the most suited framework to calculate thermodynamical properties as well as 
correlation functions. Unfortunately, Bethe ansatz is still rather ineffective to provide 
the entanglement entropies both in the coordinate [14] and algebraic [15] approaches. 
Also the special combinatoric features at A = 1/2 only allowed the calculation of 5„ 
up to ^ = 6 [TB]. By contrast, the CTM solution is very simple and the "single-particle 
levels" of the resulting i^cTM are exactly known [T3] 

Cj = 2je with e = arccoshA . (10) 

The correlation length of the model is also exactly known as function of A [TU] , but for 
what follows we are only interested in its universal part in the critical regime ^ ^ 1 

Ine^ 1^+0(6°). (11) 

The Renyi entropies are written straightforwardl}{§] 

1 



"Jn — 



oo 



J2 log(l + e-^''n - n^ log(l + e-^'') 



=0 



(12) 



This formula is exact, but does not directly allow an expansion in terms of powers of 
the correlation length, i.e. for small e. One could be tempted to perform an Euler 
Mac-Laurin summation, that gives the correct leading behavior |4j. Unfortunately 
the subleading corrections we are interested in are of the form ^^" ex e""'^ /^^, that 

§ Eq. ( |12[ | has already been derived (actually the double of it) for the XY spin chain by Pranchini 
et al. in [17) by using previous results |18) based on integrable Predholm operators and Riemann- 
Hilbert problem. It turns out that in the XXZ model, only the dependence of e on the Hamiltonian 
parameters is different, but Eq. | |12^ is identical. This could not have been known a-priori. Deriving 
this result with CTM methods allows to generalize the treatment to other models (as the XXZ spin- 
chain) for which the Riemann-Hilbert problem is not known. In Ref. |17| S„ has been rewritten in 
terms of elliptic functions, while here we use a different strategy. 
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are not obtainable by the asymptotic expansion in power of e^^ given by the Euler 
Mac-Laurin formula. One can however check straightforwardly that all the analytic 
corrections in e indeed vanish (this has been observed in [inj for n — )• oo). 

We obtain the truly asymptotic expansion for small e by using the Poisson 
resummation formula (as similarly done for Si in |20p. This formula tells us that 



where 



oo 

E 

i=-oo 



f{y) 



fim 



E/ 



2nk 
e 



f{x) coa{yx)dx . 



For Eq. (12), denoting with 

fn{x) = l0g(l 

we can rewrite the sum as 



b„ 



1 



J2ifn{ej)-nh{ej)) 



3=0 

where we used /„(0) = 
/n(y) 
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In 2. The cosine-Fourier transform of /„(x) is 
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(13) 



(14) 



(15) 



(16) 



(17) 



Plugging the last equation in the Poisson resummation formula we have that the two 
terms n/y'^ simplify, leading to the compact result 
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The term fc = gives the leading diverging term in the limit e 
can be rewritten as 



In 2 



b„ 
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(18) 
0, and so the sum 

(19) 



where the leftover sum over k is vanishing in the limit e — >■ 0. This equation is exact 
for any e as it is Eq. ( 12 ) from which we started and it is just another expansion 
of the elliptic function found in Ref. [TT]. However, in this form 5„ is ready for the 
asymptotic expansion close to e = 0. For small e, we have 



csch- 



'-k 



en 



2exp(- 



-)• 



Using Eq. ( 11 ), we have that each of these terms gives a correction of the form 



exp(- 



-)^C 



-2k/r 



(20) 



(21) 



These subleading corrections clearly depend on the precise definition used for the 
correlation length. Thus the universal character is only in the exponents 2k/n, the 
amplitudes being normalization dependent. These corrections agree with the scaling 
form proposed in the introduction with x — 2. However, for the gapless case a 
multiplicative logarithmic correction is also introduced by a marginal operator |21) . 
but it does not appear in the gapped phase. This result agrees with the recent ones 
in Ref. [211. 
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2.2. Other models 

Several other models can be treated with the help of CTMs and the same procedure as 
above for the XXZ model applies. This generalization is made particularly simple by 
the fact that in the considered cases the only change is the expression of ej as function 
of the Hamiltonian parameters. 

The most studied example is the Ising model in a transverse magnetic field with 
Hamiltonian 

^/ = -E['^>f+i+K]- (22) 

3 

This models displays a quantum critical point at ft, = 1 separating a ferromagnetic 
phase (h < 1) from a quantum paramagnetic one (h > 1). Rcnyi entropies for this 
model have been already calculated in Ref. il7] and so the following results are not 
new, but just an asymptotic expansion of known expressions. 
The "single-particle levels" of the resulting Hqtm are [13] 

r(2j + l)e for/i>l, . K{Vl-k^) 

ej — S n ■ r 1 1 With e = tt -— — , (23) 

■^ \ 2je forft<l, K{k) ^ ' 

where K{k) is the complete elliptic integral of the first kind |23], and k = min[/i, h~^]. 
An expression for the correlation length is exactly known [lOj . but for what follows 
we are only interested in the behavior close to the critical point 

ln^=y+0(6"), (24) 

valid in both phases. 

In the ferromagnetic phase h < I, the dependence of ej on j is the same as in 



the XXZ model, only e is different. Thus Eq. (19) is valid also for the Ising model 



When plugging Eq. ( 24 1 into this expression, we get a leading term corresponding to 
the central charge c = 1/2 and corrections going like ^~'=/", in agreement with the 
presence of the energy operator of dimension a; = 1 as in the conformal case E] . 

In the paramagnetic phase ft, > 1, the calculation is equivalent with the difference 
that the generalized Poisson rcsummation formula 

/(K^J +«)!) = I E /(^)e2-W., (25) 



oo 



j — -~-oo k— — oo 



must be used. From Eq. (23), we can choose a = 1/2 and b = 1 (or any multiple by 
changing the definition of f{x)). The starting formula is 

Sn = -^ E (/"(^(^' + V2)) - nfMj + 1/2))) 

= It^ E (/n(kO- + 1/2)1) -nA(|e(j + 1/2)1)), (26) 

2 1 — n ^-^ 

j=-oo 

where we notice the absence of the additive In 2 term (connected with a non-degenerate 
ground state). Thus, after Poisson resumming we have 



oo 



fc — — oo 



— -csch^^^ — — -csch 

4fc e ik en 



(27) 
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where the main difference compared to the previous case is the factor (—1)'^. Thus in 
the expansion in terms of the correlation length ^ the same power-laws enter but with 
alternating amplitudes. 

A particular simple further generalization is the XY chain for which all the 
formulas above are the same (see e.g. [SJ) and only the expression of k in terms 
of the h and 7 changes according to 

f^/W+j^^l/j forh>l, 

Xl/y/h^+l'-l foTh<l, ^^^^ 

valid for 7^ + /i^ > 1. This is reminiscent of the fact that a.t h = 1 the transition is 
always in the Ising universality class and a simple rescaling can absorb the values of 
7. Physically more interesting is the limit h = 0, not included in the case above, but 
with CTMs known [10]. One has ej — je and fc = (1 — 7)/(l + 7). The derivation 
of Sn is straightforward and one recovers the same corrections as for the Ising model, 
with a doubled leading term reflecting the values of the central charge c = 1. This is 
connected also to an exact relation between Ising and XX models [5S]. 

Also XYZ chains (see [55] for 5*1) and spin k/2 analogue of the XXZ quantum 
spin chain (see |20j for S'l) can be simply obtained from the formulas above. 

2.3. The single-copy entanglement 

The limit n — > 00 requires a separated analysis, because in Eq. ([5| the exponent of 
the corrections goes to zero, signaling the appearance of logarithmic corrections to the 
scaling as for the conformal case, ^oo is called single copy entanglement [27] [551 HH] 
and it is obtained by taking the limit n — )■ 00 before the limit e — >■ 0, since the two do 
not commute. It is easy to get the logarithmic corrections exactly. For n — >■ 00, Eq. 



(19 1 becomes 



J^ 



In 2 



00 



24e 2 -^ 

fc=i 



e 1 n^k 

2n^ " 2k ~ 



(29) 



00 



TT^ n2 e ^ 1 , Ti'k , , 

= + + > — csch . 30) 

24e 2 12 -^ 2fc e ^ ' 

fe=i 

The left over sum gives standard power-law corrections and they do not involve 
'unusual' exponents. The unusual part is in the logarithmic corrections encoded in 
the term e/12. Indeed using In^ = 7r^/(ae) (with a — \ for Ising ferromagnetic and 
a = 2 for XXZ), we have 

S^-S'Z^^ 12a log r ^^^^ 

Notice that when written in terms of e this expression has a unique correction, that 

is particularly simple for the XXZ chain since e = arccoshA. 

For the Ising model in the paramagnetic phase /i > 1, the terms (—1)*^ change 

the 1/e term, and the final result is 

7r2 e ^ (-1)'^- , ^2fc 

5*00 = :-—- — - > ^—r^csch . (32) 

24e 24 ^^ 2fc e ^ ' 

Notice that the amplitudes of the corrections are different in the two phases and their 
numerical values are of the same form of what found for the critical XX spin-chain 
[^. These results for the single copy entanglement in XY chains have been reported 
also in [HIIH]. 
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3. Relation to Virasoro characters 

The above examples may be considered as special cases of a more general phenomenon: 
that in many integrable models satisfying the Yang-Baxter relations, traces of 
powers of the CTM may be expressed in terms of the Virasoro characters of the 
corresponding CFT of the critical theory. These enjoy simple properties under modular 
transformations which allow us to extract the 'unusual' corrections to scaling discussed 
earlier. 

To be specific, the Boltzmann weights of many integrable models lie on a curve 
specified by an elliptic modulus q, where the critical point corresponds to q -^ 1. The 
examples of quantum spin chains considered above correspond to q = e~^'^. Baxter 
[TU] showed that in these models the eigenvalues of the 4th power A^ of the CTM are 
all of the form q"-^+b ^ for integer iV > 0, with a computable degeneracy d^- It was 
then observed in the late 1980s (see, for example, Refs. [291 [30]) that the degeneracy 
factors are just those appearing in the highest weight representations of the Virasoro 
algebras in the CFT which describes the scaling limit at the critical point, that is 

XA(<?) = g-^/^^+^f:d^g^, (33) 

where c is the central charge and A is the highest weight of the representation, giving 
the dimension x = 2A of a scaling operator of the theory. The value of A depends 
on the value of the spin at the origin, and the boundary conditions at infinity. This 
observation, although not proven in general, was based on the facts that (a) in the 
scaling limit at the critical point, but in a finite annulus with outer and inner radii i?> 
and i?<, the eigenvalues of the CTM have precisely this form, with conjugate modulus 
q (defined below) equal to R^/Ry - this follows from boundary CFT [3T]; and (b) in 
the integrable RSOS models, the eigenstates of the CTM correspond to a semi-infinite 
walk on the corresponding Am diagram, and the Virasoro generators act simply on 
the space of these paths. (In fact it was shown by the Kyoto group [35] that for higher 
integrable models obtained by fusion the result is not always a simple character but 
rather a branching coefficient of an affine algebra. However these still enjoy simple 
modular properties as described below.) 

This means, in particular, that the partition function TrA^ is proportional to 
XA{q'^) for some A (which depends on, for example, which massive phase the model is 
in). Without loss of generality we take a = 1. In general the partition function may be 
a linear combination of characters, depending on the choice of boundary conditions at 
the origin and infinity. This complicates the argument without changing the general 
conclusion, and we shall assume that the boundary conditions are such as to pick out a 
single character, as in the examples discussed earlier. However we remark that, in the 
context of entanglement entropy, we should specify if we fix the spin at the origin in 
the direct or in the dual lattice. In all examples above, we fixed the dual variable, by 
leaving the bond between the two halves free. However, while at first it could sound 
strange, also fixing an actual spin at the origin makes sense for the entanglement 
entropy, because in the CTM this does not correspond to fixing it in the hamiltonian 
(which would effectively divide the chain into two non- interacting halves). Rather it 
means projecting the ground state |0) into a subspace in which the spin at the origin 
is fixed, and measuring the entanglement between the two halves of the chain in this 
subspace. Thus we expect to be able to explore other values of A by such a procedure, 
and will therefore keep it general in what follows. 
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The point now is that the trace of A*" is, in the basis in which the CTM is 
diagonal, given simply by replacing q -^ q", and thus 

Trp'A^r^T^- (34) 

{XA{q)) 

We are interested in the limit q — > 1 close to the critical point, where the series 



expressions ( 33 ) for XA (q) are not very useful. However |33j , the Virasoro characters 



transform linearly under a modular transformation g — > q, where, if q =^ e , 

XA (g) ^^S^ XA' (q) , (35) 

A' 

where S^ is the modular S-matrix which characterizes the CFT. Notice that the 
modular invariance plays exactly the same rule as Poisson resummation above, but it 
is more general. 

As 5 —> 1, (7 —> 0, so it is straightforward to extract the leading behavior and all 
the corrections. The leading term comes from A' = 0, in which case we see that 



TTp^X-is^AY (q 



\ 1-n f^-c/24\ " 



which gives 



Sr,^~^(l + ^]\nq + In (5^) + o(l) , 



(36) 



Comparing with Q, we see that <? oc ^ ^. The 0(1) term is related to the AfHeck- 

Ludwig boundary entropy 34 g^ = In (S'a) • Note that we get twice the boundary 
entropy because there are two semi-infinite chains adjoining the spin at the origin. The 
examples discussed earlier presumably correspond to A = and the corresponding 
0(1) term can be absorbed into the non-universal constant of proportionality between 
q and ^~^. 

The corrections now come from (a) integer powers of q, that is powers of ^^^, in 
the expansion of the character, and (b) other characters with A' > 0. The leading 



terms of the latter form then give rise to corrections to ( 36 ) 



(i--r'E#(9^'^"--9^^') + ---- (37) 

A'>0 '^ 

(The second power is less important for n > 1 but is included so as to ensure the 
finiteness as n — > 1.) Recalling that 2 A' = x' , we see that we get corrections of the 
form ^~^ ' " as expected. However, unless the matrix element S^ happens to vanish, 
in principle we get such contributions from all the scaling dimensions x' of primary 
fields, not just the relevant ones. Also, the higher order corrections in general involve 
powers which are all possible integer linear combinations of x' /n and x'. The evidence 
of these further corrections has been reported for simple gapless spin chains [8] , but 
it would be as interesting to report them in the gapped case. 

We can also derive a general result for the single-copy entanglement Sao- Using 



Eq. ( 34 1 we see that 

5oo = -(-c/24 + A)ln(7 + lnxA(g). (38) 
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The second term, after using the modular transformation, goes Uke (— c/24 + A') \ogq 
plus power-law corrections. We expect the leading term in the sum over A' to come 
from A' = 0. Thus we get a universal term (c/12)ln^. This corresponds to the first 
term in Eqs. (pO 
discussed in Sec. 



) an d (32). The first term in Eq. (381 gives the 'unusual' correction 
2.3 Since q ex ^^^, \nq ^ — 27r^/ln^, so this correction has the form 



^-^(-f^/ + ^) . (39) 

This agrees with the Ising case in the ferromagnetic phase given by Eq. (|31[) using 



c — 1/2, A = 1/16, and in the paramagnetic Eq. (32) using c — 1/2 and A = 



4. Concluding remarks 

We showed that the corrections to the scaling of the Renyi entanglement entropies 
in gapped systems display the universal form ([5| in the case when an infinite line is 
divided into two semi-infinite subsystems. We provided few explicit examples using 
known corner transfer matrices and we argue that in the general case these corrections 
are a consequence of the modular invariance of the traces of powers of the CTM when 
expressed in terms of the Virasoro characters. 

This result generalizes straightforwardly to more complicated bipartitions 
provided the correlation length ^ is smaller than all separations. Eq. (Isl) only gets 
multiplied by the number of boundary points between A and B. The rigorous result for 
a finite interval of length £ obtained for the XY model ^7] confirms this. However, for a 
finite interval, only when £ becomes much larger than ^, the asymptotic behavior would 
be visible. For smaller values of £, a complicated crossover between the conformal 
result and the asymptotic one takes place as already known for the leading terms |35) . 
The characterization of all these corrections is not only an academic task. In the 
case of the entanglement of two disjoint intervals in a conformal model, their precise 
knowledge has been fundamental to recover the CFT predictions ^E[ [37] in numerical 
calculations [38l |39l HQI HI]. We believe that the same is true for massive systems, 
especially in the case of non-integrable models when numerical calculations are the 
only way to attack the problem. A final question is whether these unusual corrections 
are also present for other entanglement measures that displayed them in the gapless 
phases as the valence bond entanglement [12], and if yes whether it is possible to 
calculate them exactly. 
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